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We show that spin-orbit coupling (SOC) gives rise to pairing instability in a highly polarized 
two-dimensional Fermi gas for arbitrary interaction strength. The pairing instability can lead to a 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)-like molecular state, which undergoes a first-order tran- 
sition into a pairing state with zero center-of-mass momentum as the parameters are tuned. These 
pairing states are metastable against a polaron state dressed by particle-hole fluctuations for small 
SOC. At large SOC, a polaron-molecule transition exists, which suggests a phase transition between 
the topological superfluid state and the normal state for a highly polarized Fermi gas in the ther- 
modynamic limit. As polarization in a Fermi gas with SOC is induced by the effective Zeeman field, 
we also discuss the influences of the effective Zeeman field on the ground state of the system. Our 
findings may be tested directly in future experiments. 
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Spin-orbit coupling (SOC), a non-Abelian gauge field, 
has been shown to play a fundamentally important role 
in many interesting systems in condensed matter physics, 
e.g., topological insulators [1], quantum spin Hall materi- 
als [2], etc. The recent development of a synthetic gauge 
field, SOC in particular, in ultracold atoms has stimu- 
lated a tremendous amount of interest in the study of 
the effects of SOC within these systems [3-5]. By break- 
ing the inversion symmetry, the SOC may induce novel 
quantum phases, e.g., the unconventional superfluidity 
in an ultracold Bose gas [6, 7], or the topological super- 
fluid phase (TSF) in a polarized Fermi gas with s-wave 
pairing order [1, 8, 9, 11-13]. Notably in the latter case, 
Majorana zero modes can be stabilized at the center of 
vortex excitations, which may be used as resources for 
topological quantum computation [14, 15]. 

For a polarized two-dimensional (2D) Fermi gas with 
SOC, a peculiar behavior is the existence of the pairing 
instability in the large polarization limit [1-3] . This is in 
stark contrast to the case without SOC, where the pairing 
state becomes unstable against a normal gas beyond the 
so-called Chandrasekhar-Clogston limit [18]. The per- 
sistence of the pairing instability can be attributed to 
the breaking of inversion symmetry, which modifies the 
topology of the Fermi surface and renders singlet s-wave 
pairing possible in the large polarization limit [1]. This 
can lead to the interesting scenario where pairing exists 
in the presence of a single Fermi surface, i.e., when the 
chemical potential lies in the gap between the two helic- 
ity bands. In this case the system is in the TSF phase. 
For a uniform 2D gas, the TSF phase is stable in the 
large polarization limit on the mean-field phase diagram 
[19]. However, it has been shown theoretically [5, 20] and 
demonstrated very recently in experiments [22], that in 
the absence of SOC the ground state of a highly polar- 



ized 2D Fermi gas is a polaron state in the weak-coupling 
limit, i.e., an impurity atom dressed by particle-hole fluc- 
tuations of the Fermi sea. In the presence of SOC, nat- 
urally one expects the interplay of SOC, pairing, and 
fluctuation leads to rich physics. In particular, it is in- 
teresting to study the stability of the TSF phase against 
a normal state with particle-hole fluctuations in a highly 
polarized Fermi gas. 

In this work, we investigate a spin-orbit coupled 2D 
Fermi gas in the large polarization limit. In the presence 
of SOC, the spin polarization can be induced by an effec- 
tive Zeeman field, which is tuned by adjusting the laser 
parameters in a typical scheme for synthetic SOC [3, 4]. 
To model the ground state of the system, we adopt varia- 
tional ansatz states following Refs. [5, 20, 23, 24], which 
effectively project the wave functions into the subspace of 
the large polarization limit. We then study in detail the 
properties of both the molecular state and the polaron 
state in the presence of an effective Zeeman field. We 
find that under appropriate effective Zeeman fields the 
SOC-induced pairing instability leads to a Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO)-like pairing state [25] with 
nonzero center-of-mass momentum in the weak coupling 
limit. The FFLO-like pairing state can undergo a first- 
order transition into a pairing state with zero center-of- 
mass momentum as SOC increases or as the interaction 
is tuned. For small SOC, we find that the energy of a 
polaron state is always lower than that of the molecu- 
lar pairing state. However, a polaron-molecule transition 
exists for sufficiently large SOC. This suggests a phase 
transition between the normal state and the TSF state 
in the thermodynamic limit. Finally, we show that the 
boundaries between the different states can be shifted as 
the effective Zeeman field increases. With progress in 
the experimental investigation of 2D Fermi gases [22, 26] 
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and the recent realization of SOC in a degenerate Fermi 
gas [27] , our study has interesting implications for future 
experiments. 

Model.- We consider the system in the large polariza- 
tion limit induced by an effective Zeeman field h. In the 
presence of Rashba SOC, the Hamiltonian can be written 
as 
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where e k = h 2 k 2 /2m, the pseudospin of the atoms a = (f 
a is the SOC strength, </? k = arg (k x + ik y ), and S is 
the quantization volume in two dimensions. Here we only 
consider the case where atoms of different spin species 
have the same mass m. The bare s-wave interaction rate 
U should be renormalized following the standard relation 
in two dimensions [28] , 
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where E b is the binding energy of the two-body bound 
state in two dimensions in the absence of SOC. For 
cold atom systems, this two-body binding energy can 
be tuned, for instance, via the Feshbach resonance tech- 
nique. 

Molecular state - We first investigate the pairing insta- 
bility in the presence of SOC, and the properties of the 
resulting molecular state. Consider a variational ansatz 
of molecular state with center-of-mass momentum Q: 
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where |iV — 1) represents a Fermi sea with N — 1 spin- 
up atoms. Due to the SOC, we now have triplet-pairing 
components |<^, m addition to the singlet-pairing 

wave function </>^. The momentum subscripts in the 
wave functions are constrained by Pauli blocking, such 
that k > kp for <j>^ and , and |Q — k| > kp 
for <fiU , where kp is the Fermi wave vector given by 
Ep = h 2 k F /2m, and Ep is the Fermi energy of spin-up 
fermions. To focus on the properties of the pairing insta- 
bility, we neglect the particle-hole fluctuations here and 
only consider the 'bare' molecular state in Eq. (3). Note 
that by taking the ansatz in Eq. (3), we have effectively 
projected the ground state into a subspace that corre- 
sponds to the large polarization limit where few spin- 
down atoms coexist with a polarized Fermi sea of spin-up 



atoms. Terms with further spin flips are suppressed by 
the effective Zeeman field h and are projected out. Here, 
to be of more experimental relevance, we fix the effec- 
tive Zeeman field, which typically corresponds to fixing 
the laser parameters for synthetic SOC. As a result, the 
population of the spin-down atoms in the ground state 
fluctuates slightly around unity, given appropriately cho- 
sen Zeeman field strengths. 

Minimizing the functional (M\ (H — Em) \M), we get 
a self-consistent equation for the ground state energy Em 
of the molecular state, relative to the Fermi sea of TV 
spin- up atoms [29]: 
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where 8(x) is the Heaviside step function. To find the 
molecular ground state, we should further minimize the 
energy solved from Eq. (17) with respect to the center-of- 
mass momentum Q. In the absence of SOC, the center- 
of-mass momentum of the molecular state is zero for large 
Eb, becomes finite, i.e. FFLO-like, at E b = 2Ep through 
a second-order transition, and approaches the Fermi wave 
vector kp at E b — 0.5E F where the bound state merges 
into the continuum [20] . Importantly, there is no pairing 
instability for Ej, < 0.5Ep without SOC in the large po- 
larization limit. We will see that this simple picture is 
drastically modified by SOC. 

We numerically solve Eq. (17) and minimize the so- 
lution with respect to the center-of-mass momentum Q 
for a fixed Zeeman field h = 0.5Ep . The evolutions 
of the molecular energy as well as the center-of-mass 
momentum Q are shown in Fig. 1. Compared to the 
case without SOC, an outstanding difference is that the 
pairing instability persists into the weak-coupling limit 
even with infinitesimally small SOC. The bound state en- 
ergy approaches an SOC-dependent asymptotic value in 
the weak-coupling limit and never crosses the continuum 
threshold Em — h — a 2 k 2 F /2Ep. This is consistent with 
the previous many-body calculations, where the super- 
fluid order parameter does not vanish for any polariza- 
tion in the presence of SOC [1-3]. The pairing instability 
in the large polarization limit is actually a consequence 
of the combined effects of the SOC and the effective Zee- 
man field. While SOC breaks the inversion symmetry 
and mixes the spins into different helicity bases, the ef- 
fective Zeeman field breaks the time reversal symmetry 
and opens a gap in the energy spectrum. As the spins 
are mixed in the lower branch, s-wave pairing is possi- 
ble even in the large polarization limit [1]. Numerically, 
the pairing instability shows up as a singularity in Eq. 
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FIG. 1: (Color online) (a) The molecular energy (Em — 
h)/EF as a function of Eb/EF for different SOC strengths 
and a fixed effective Zeeman field h/EF = 0.5: akp/EF = 
(thin solid line), akp/EF = 0.2 (thin dashed line), akpEF = 
0.6 (bold solid line), akp/EF = 1 (bold dashed line), (b) 
Evolution of the center-of-mass momentum of the molecular 
ground state for different SOC strengths and a fixed effective 
Zeeman field h/Ep — 0.5. 



(17) when Eb — > 0, similar to the behavior of the gap 
equation in the many-body case. From the many-body 
mean-field calculations, the ground state of a uniform 2D 
Fermi gas is a TSF phase in the large polarization limit 
[19]. Apparently, the pairing state that we discuss here is 
related to the topological superfluid phase in the thermo- 
dynamic limit when a finite density of spin-down atoms 
are present [2]. 

Another important observation is that, for small SOC, 
the pairing instability leads to an FFLO-likc pairing state 
with a finite center-of-mass momentum Q whose magni- 
tude approaches the Fermi wave vector kp in the weak- 
coupling limit. This is shown in Fig. 1(b). As the 
interaction is tuned toward the strongly interacting re- 
gion, the center-of-mass momentum becomes smaller and 
drops to zero at a critical E^, where the system under- 
goes a first-order transition to a pairing state with zero 
center-of-mass momentum. The critical E^ for this first- 
order transition is, in general, a function of the SOC 
strength and the effective Zeeman field h. As demon- 
strated in Fig. 2(a), for small h, E£ decreases monoton- 
ically as the SOC strength becomes larger, and eventu- 
ally vanishes at a critical SOC strength a c . Beyond a c 
only pairing states with zero center-of-mass momentum 
are stable for arbitrary interaction strength. For larger h, 
the dependence of E% on the SOC strength becomes non- 
monotonic. This suggests that the pairing state becomes 
non-FFLO-like in the weak coupling limit for large h, and 
that one may encounter up to two first-order transitions 
if the interaction is tuned with fixed SOC and Zeeman 
field. Hence for a fixed Zeeman field, the stability of the 
FFLO state is affected by SOC in two contrasting ways: 
while the competition between the FFLO state and the 
pairing state with zero center-of-mass momentum limits 
the stability region of the FFLO state, the pairing insta- 
bility induced by SOC effectively increases the stability 
region of the FFLO state, especially in the weak-coupling 
regime. An outstanding question here is what are the 




FIG. 2: (Color online) (a) The critical El for the first-order 
transition to FFLO-like pairing state as a function of the SOC 
strength for various effective Zeeman fields: K/Ef — 0.5 (solid 
line), H/Ef — 0.95 (dashed line), h/Ep — 1.1 (dash-dotted 
lines), (b) The critical SOC strength in the weak-coupling 
limit Eb = as a function of the effective Zeeman field h. 

topological properties of the FFLO pairing state, which 
deserves further investigation in the future. 

Now that we understand the nature of the pairing 
states in the weak-coupling limit, we may derive an ana- 
lytical expression for the molecular energy in this limit. 
Similar to the divergence of the gap equation in the 
many-body case [1, 2], the summation in Eq. (17) must 
diverge when Eb — > 0, as dictated by the pairing insta- 
bility. The energy of the molecular ground state in the 
weak-coupling limit corresponds to the lowest-lying sin- 
gularity of the summand on the right-hand side of Eq. 
(17). To evaluate this energy, we set the denominator in 
the summand to zero and minimize the molecular energy 
with the constraint k > hp. As the center-of-mass mo- 
mentum Q is either kp or in the weak-coupling limit, we 
need only consider these two cases. For the FFLO state 
with Q — kp, we find the lowest energy for the occur- 
rence of the sing ularity at E S M = (3h- yjh 2 + 8a 2 k 2 F )/2 1 
with {(p — >0,k — > kp} (cp is the angle between Q and 
k). For Q = 0, the asymptotic energy is given by 
E f M = h-E F -(h 2 +a 4 k F /E F )Ep/2a 2 k F for a l k%/E\- 
h 2 > Aa 2 k 2 Fl and E S M = E F + h - ^Jh 2 + Aa 2 k 2 F for 
a 4 k%/E 2 F - h 2 < 4a 2 k 2 F . With these, the critical SOC 
strength a c at which the FFLO-like pairing state vanishes 
at Eb = can be calculated as a function of the effective 
Zeeman field h [see Fig. 2(b)]. Apparently, the stabil- 
ity region of the FFLO-like pairing state decreases as the 
Zeeman field h increases. In particular, the FFLO-likc 
pairing state vanishes in the weak coupling limit beyond 
h/Ep = 1 for arbitrary SOC . 

Polaron state.- The molecular states that we have dis- 
cussed are the ground state of the system if only the 'bare' 
polaron state \P) — i/j Cq^\N) is considered [5]. In two 
dimensions, we should include the particle-hole fluctua- 
tions above this 'bare' polaron state for a more realistic 
calculation. The energy of the dressed polaron state is 
then lowered and the molecular state in Eq. (3) may no 
longer be the ground state. Indeed, in the absence of 
SOC, it has been shown previously that for a single im- 
purity atom in the presence of a polarized Fermi sea in 
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FIG. 3: (Color online) (a) The energy of the polaron state 
as a function of Eb/Ep for various SOC strengths with fixed 
effective Zeeman field h/Ep = 0.5: akp/Ep = (thin solid 
line), akp/Ep = 0.2 (thin dashed line), akpEp = 0.6 (bold 
solid line), akp/Ep — 1 (bold dashed line), (b) The polaron- 
molecule transition boundary on the {Eb/Ep, akp/Ep} plane 
for various Zeeman fields: h/Ep = 0.5 (solid line), h/Ep = 
0.95 (dashed line), h/Ep = 1.1 (dash-dotted line). 

two dimensions, a polaron state dressed by particle-hole 
fluctuations has lower energy than the 'bare' molecular 
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where E^ q = e q _k + e k — e q , and the third term on the 
left-hand side corresponds to an energy shift due to the 
SOC-induced particle-hole fluctuations in the Fermi sea 
of spin-up atoms. In the weak-coupling limit, this should 
be the dominant contribution as the summation on the 
right-hand side of Eq. (6) vanishes at Eh = 0. Hence 
the polaron energy in the weak-coupling limit is E p = 
(3ft.— \Jh 2 + 4-7ra 2 fc|,) /2. Comparing this energy with the 
molecular energy in the same limit, we can solve for the 
critical SOC strength of the polaron-molecule transition. 

We numerically evaluate the energy of the polaron 
state, and calculate the critical SOC for the polaron- 
molecule transition for various h as the interaction is 
tuned. As is clear from Fig. 3, for sufficiently large 
SOC strength, molecular states are always favored, de- 
spite the shift in polaron energy by SOC. In addition, 
we find that the critical SOC for the polaron-molecule 
transition is always greater than the critical SOC for the 
FFLO-like pairing states under the same effective Zee- 
man field. Therefore at the polaron-molecule transition, 
the center-of-mass momentum of the molecular state is 
always Q = 0. As the molecular state with zero center- 
of-mass momentum corresponds to the TSF state in the 
mean field diagram, this polaron-molecule transition im- 
plies a phase transition between the TSF state and the 
normal state for a highly polarized 2D Fermi gas in the 



state that is not dressed by such fluctuations [5, 20]. 

We consider a variational ansatz for such a polaron 
state with one particle-hole density fluctuation in the 
large polarization limit 

\P) = ik^N) + £ V k t q e q _ k4 c k)tCq , t |A0 

k>k F ,q<k F 

+ E ^A-^A,^M N ) 

\t,q<k F 
q<k F ,k>k F 

(5) 

Similar to the molecular case, the effective Zeeman field 
is applied to induce the spin polarization. Again, we fix 
the Zeeman field in the following for simplicity and for 
experimental relevance. 

Minimizing (P\ (H — Ep) \P), we find a self-consistent 
equation for Ep, the energy of the polaron state 
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I 

thermodynamic limit. 

For a 2D Fermi gas, fluctuations are important, and a 
'bare' molecular state in Eq. (3) may not be accurate. 
It has been pointed out in Ref. [5] that if particle-hole 
fluctuations are included in the molecular states as has 
been done for the polaron state, the energy of the molec- 
ular state will be lowered. In the presence of SOC, we 
expect a similar scenario will shift the critical point of 
polaron-molecule transition to smaller SOC strengths. 

Conclusion- We have studied the novel physics in- 
duced by SOC in a 2D Fermi gas in the large polarization 
limit generated by an effective Zeeman field. With the 
interplay between SOC, pairing, and particle-hole fluc- 
tuations, the system exhibits many interesting proper- 
ties, e.g. SOC-induced pairing instability, FFLO pairing, 
polaron-molecule transition, etc. In particular, our re- 
sults suggest that the particle-hole fluctuations have con- 
siderable impact on the stability of the TSF state and can 
modify the mean-field phase diagram in the large polar- 
ization limit. As polaron-molecule transitions in a highly 
polarized Fermi gas have been experimentally probed re- 
cently both in three dimensions [30] and in two dimen- 
sions [22], our results may be tested directly in future 
experiments. 
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Supplementary material 



In this supplementary material, we present in detail the variational formalism adopted in the main text. We wish 
to study the system in the large polarization limit, where a pairing instability persists in the many-body calculations 
[1-3]. In the presence of SOC, the spins are mixed, and the large polarization limit can be achieved by imposing an 
effective Zeeman field which prevents spin- flipping. The implementation of this effective Zeeman field depends on the 
specific scheme in realizing the synthetic SOC. For the scheme in Ref. [4] for example, the effective Zeeman field is 
tunable by adjusting the Rabi-frequency of the Raman lasers. In the canonical ensemble with a total of N particles, 
the large polarization limit corresponds to the subspace with few spin-down atoms in the presence of a fully-polarized 
Fermi sea of spin-up atoms. As systems with different population of spin-down atoms can be connected to one 
another by a simple scaling of the parameters, we further specify the number constraint in the large polarization limit 
so that there is one spin-down atom on average. To model the ground state in this large polarization limit, we adopt 
variational ansatz states which effectively project the wave functions into this subspace. As we will see in the following, 
terms with more than two spin-down atoms are left out of the ansatz, as their contributions are typically small due 
to the presence of the effective Zeeman field. In practice, the Zeeman field h can be self-consistently determined by 
imposing the number constraint. Here, instead of fixing the average number of spin-down atoms, we fix the effective 
Zeeman field h, which is of more experimental relevance. The average number of spin-down atoms in our calculation 
therefore fluctuates around unity. 



Molecular state 

In this section, we present the formalism for the molecular state. Consider an interacting Fermi gas with SOC and 
an effective Zeeman field h. The Hamiltonian is 

H = E £ k4,a C k^ + U E c k,T C k',; c k'+q,iCk-q/t + E ( akc lt C ^ + a * fcc L c k,t) + ^E a k,4 a k,4> ( 7 ) 
k,<r k.k' ,q k k 

with ek = h 2 k 2 /2m. The bare interaction rate U should be renormalized following the standard recipe in 2D 

U = ~S^ E b + 2e^ (8) 

k 

where Ej, is the two-body bound state energy in 2D in the absence of SOC. 
For the molecular state, the variational wave function can be written as 
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Q = E (^ c Q-k^k,t + ^ c Q-k,t c k,T + ^ c Q-k,+ c L) 1^-1). ( 9 ) 



where \N — 1) represents the Fermi sea with N — 1 spin- up fermions. Importantly, due to the presence of the Fermi 
sea, the coefficients are constrained such that k > Uf for </>j^ and </>j^; and |Q — k| > kp for <^>£ . The normalization 
condition for the total wave function is 

<*m|*m) = E i^i 2 + E (iQ- k i - m (i4 T i 2 - tfVff-k) + E (i^i 2 - 

k>kjp k>kf k 

= i, (io) 

where 0(x) is the Heaviside step function. Note that the last two terms can be recast into a more symmetric form 

E O^l 2 - <^Q-k) = E \ " ^Q-k) Wk - 0Q-k) • (11) 



We may then evaluate the quantity ^ m\H\^ m) and minimize it to derive the equations for the coefficients 0k in the 
wave function. 
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The expectation value of the energy becomes 
(* M \H\9 M ) = E |^| 2 (e k + e Q _ k + h) + £ \ (f^ - #£ k ) - ^_ k ) (e Q _ k + e k + 2h) 0(|Q - k| - M 



k>k 



E 0(|Q-ki|- M (< T - 0Q n ki ) |Q - k t | + h. 



c. 



,k± >/cj? 



M>ki 



(12) 



The extremal conditions of the energy functional above then yield the following set of equations for the variational 
coefficients 



(e k + e Q -k + h) 4£ + U E 4' + «* IQ ~ k l #T + = ^ 



k'>k F 
1 



(ck + eq-k)^r + Qt|Q-k|^ 



1 



(e k + eq- k + 2h) <j>£ + a*fc</£ 



where we have defined the properly symmetrized coefficients 

</> k = k - ^Q- k - 



(13) 

(14) 
(15) 

(16) 



It is important to keep in mind that in the equations above, we have the constraints k > kp and |Q — k| > kp. 
Following the standard procedure, we derive the equations for the ground state energy 



l -= E 



k>ki 



1 



E - (e k + £Q- k + h) - 



2|q| 2 fc 2 



-E-(e k +eQ-k+2h) 



-0(|Q-k|-fc F )^ 



2|a| 2 |Q- k | 2 



-B-(ek+EQ-k) _ 



(17) 



Finally, the coefficients can be solved once the ground state energy is fixed 



E-B 



E 



fe>feF (E-B) 2 



where 



B = (ek + eq-k) + h 



2\a\ 2 k 2 



C = 1 + 9(\Q - k| - k F ) 



£ - 2/i - (e k + e Q _ k ) 
2M 2 |q~k| 2 



|Q-k|-/c F ) 
2\a\ 2 k 2 



2|a| 2 |Q-k| 2 

E - (e k + e Q _ k ) 



(£ - eq-k - e k ) 2 (E - e Q _ k - e k - 2/i) 2 



(18) 

(19) 
(20) 



The average number of spin-down atoms can be calculated from N± — Efe>fe F I^I^I 2 + E k 1^ 



kl ^| 2 . In Fig. 4, we 



illustrate N± as a function of E), for fixed SOC strength a and various Zeeman fields. For a fixed Zeeman field with 
appropriate magnitude, the population of the spin-down atom fluctuates around one, as expected. 



Polaron state 



Similar to that of the molecular state, the variational wave function for the polaron state can be written as 

\p)= 1 p ci l \N}+ e < q c Q-k 4 <t c «i.ti JV >+ E ^kA-kA^mn) 

k>kp ,q<kF \i.,q<kp 

+ e "(ii- k i - ^)^! q <-k,t c 3c, t ^,ti^)- 



(21) 



q<kp .k>kp- 



1.4 



1.2 




1 1.5 2 

Eb/Ep 

FIG. 4: (Color online) Average number of spin-down atoms as a function of Eb with fixed SOC strength akp/EF = 0.2 for 
various effective Zeeman fields h: H/Ef = (dashed), H/Ef = 0.5 (solid), H/Ef = 1-1 (dash-dotted). The abrupt changes 
indicate first-order transition between the FFLO-like pairing states and pairing states with zero center-of-mass momentum. 



The normalization condition for the coefficients are 



{p\p) = \M 2 + E r 

k~>kp ,q<k 



n 

k.q 



E (K 



E edq-ki-^n^l -0£k, q 

q<kp ,k>kp 

: 1. 



(22) 



Minimizing the energy functional in the presence of the Zeeman field, we have 



Uj2^o + U E ^& + E a ^qtq = E ^o - hiPo 

q<kjp k>kp ,q<kp q<kp 

Ek^H + U ^ + u E ^', q + a ^k! q + «* k - ki = B^| q - 



k'>k f 



n 



1 " l -'tk, q = ^^k, q . 



(23) 
(24) 

(25) 
(26) 



where i? fcg = e k + e q -k - e q , and 



— <Pk, q — <P q -k, q - 



(27) 



The equation for the energy is 



q<k B 



2\a\ 2 q 2 
E-2h 



= E 

q<k F 



1 

E7 



E 

k>k F 



E — Ek q — h — 



\a\ 2 k 2 
±(E-E kq -2h) 



-0(|q-k| 



(28) 



In the absence of a and this equation reduces to the corresponding equation for the polaron energy in Ref. [5]. 
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From the renormalization condition and the eigen energy equations, we have 



where we have defined 



1+ J2 c ^ 

k>kjp ,q<kp 

B, 



E-A 



kq 



2ak 



n 

kq E - E kq - 2h ^ q 

n 



2a|q-k| n 

Vk,q> 



k,q E — E, 



Akq = Ekq 



B n = 



kq 



2a 2 k 2 



E — E kq — 2h 



0(|q-k|-fc F ) 



2a 2 |q-k| : 



feg 



E7 Sfc'; 



c fc9 - 1 + 



2a 2 fc 



21,2 



{E - E kq - 2hf 



|q- k| - hi 



2a 2 |q-k| : 

(E - E kq f 



(29) 



We have checked numerically that for effective Zeeman fields on the order of Ep, the population of the spin-down 
atoms remains close to unity. 



* Electronic address: wyiz@ustc.edu.cn 

* Electronic address: wzhangl@ruc.edu.cn 

[1] S. Tewari, T. D. Stanescu, J. D. Sau, and S. Das Sarma, New J. Phys. 13, 065004 (2011). 

[2] J. Zhou, W. Zhang, W. Yi, Phys. Rev. A 84, 063603 (2011). 

[3] J.-N. Zhang, Y.-H. Chan, and L.-M. Duan, arXiv: 11 10.2241. 

[4] Y.-J. Lin, K. Jimenez-Garci'a, and I. B. Spielman, Nature (London) 471, 83 (2011) 

5 M. M. Parish, Phys. Rev. A 83, 051603 (R) (2011). 

I 



